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We make a detailed study of boson star configurations in Jordan-Brans-Dicke theory, studying both 
equilibrium properties and stability, and considering boson stars existing at different cosmic epochs. 
We show that boson stars can be stable at any time of cosmic history and that equilibrium stars are 
denser in the past. We analyze three different proposed mass functions for boson star systems, and 
obtain results independently of the definition adopted. We study how the configurations depend on 
the value of the Jordan-Brans-Dicke coupling constant, and the properties of the stars under extreme 
values of the gravitational asymptotic constant. This last point allows us to extract conclusions about 
the stability behaviour concerning the scalar field. Finally, other dynamical variables of interest, 
like the radius, are also calculated. In this regard, it is shown that the radius corresponding to the 
maximal boson star mass remains roughly the same during cosmological evolution. 
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I. INTRODUCTION 



Although boson stars are so far entirely theo- 

retical constructs, they give rise to one of the simplest 
possible stellar environments in which to study gravita- 
tional phenomena mathematically. One can find numeri- 
cal solutions which are nonsingular and yet exhibit strong 
gravitational effects. Many of their properties bear close 
resemblance to those of neutron stars. 

Boson stars were first conceived as Klein-Gordon geons 
— systems held together by gravitational forces and com- 
posed of classical fields. They are a gravitationally-bound 
macroscopic state made up of scalar bosons. As with 
neutron stars, the pressure support which leads to their 
existence is intrinsically quantum. For neutron stars, the 
pressure support derives from the Pauli exclusion princi- 
ple, and for boson stars this is replaced by Heisenberg's 
uncertainty principle. Assuming that the quantum state 
contains sufficient particles for gravitational effects to 
be important, and that particle interactions can be ne- 
glected, an estimate of the mass is readily obtained as 
follows. For a quantum state confined into a region of 
radius R, and with units given hy h — c — 1, the boson 
momentum is p — 1/R. If the star is moderately rela- 
tivistic, p m, then R ~ 1/to. If we equate R with the 



Schwarzschild radius 2M / rrip^ (recall that G 



"PI 



), we 



find Af ~ rrip^/m. 

In practice, one assumes the existence of a classical 
scalar field with a given Lagrangian density, and adopts 
an ansatz for its time dependence which implicitly en- 
codes the Heisenberg uncertainty. This time dependence 
is of course of a form which still permits a static metric. 
With these ingredients, one then solves Einstein's equa- 
tions, something which must be done numerically. When 
no self-interaction term is present in the Lagrangian den- 
sity, the masses concur with the estimate above. How- 



ever, if self-interaction is present it is typically the domi- 
nant contributor to pressure support, and leads instead to 
masses of order rrip^/m? . If the boson mass is comparable 
to a nucleon mass, this order of magnitude is compara- 
ble to the Chandrasekhar mass, about IMq Thus, 
boson stars arise as possible candidates for non-baryonic 
dark matter, and are possibly detectable by microlensing 
experiments. 

Boson stars have been widely studied in general relativ- 
ity, where the basic model has been extended in various 
ways, such as including a C/(l) charge allowing a mix- 
ture of boson and fermion components 0, or including 
a non-minimal coupling of the boson field to gravity 
These and others works are summarized in two reviews 
and more recently in Ref. The possibility of 
direct observational detection of boson stars was studied 
recently in Ref. where it was asked whether radi- 
ating baryonic matter moving within a boson star could 
be converted into an observational signal. Unfortunately, 
any direct detection looks a long way off. 

Given the simplicity of the boson star, it is natural to 
examine boson star solutions in theories of gravity other 
than general relativity, to examine whether new phenom- 
ena arise. The most-studied class of such theories are the 
scalar-tensor theories of gravity [0, which include the 
Jordan-Brans-Dicke (JBD) theory as a special case. In 
these, Newton's gravitational constant is replaced by a 
field (p known as the Brans-Dicke field, the strength of 
whose coupling to the metric is given by a function uj{(j)). 
If w is a constant, this is the JBD theory which is 
the simplest scenario one may have in this framework. 
General relativity is attained in the limit l/w — *■ 0. To 
ensure that the weak- field limit of this theory agrees with 
current observations, lo must exceed 500 at 95% confi- 
dence ]ll| from solar system timing experiments, i.e. ex- 
periments taking place in the current cosmic time. This 
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limit is both stronger and less model-dependent than lim- 
its from nucleosynthesis . Scalar-tensor theories have 
regained popularity through inflationary scenarios based 
upon them and because a JBD model with lu = —1 
is the low-energy limit of superstring theory jist . 

The first scalar-tensor models of boson stars were stud- 
ied by Gunderson and Jensen [T^ ], who concentrated on 
JBD theory with lo = 6. This was generalized by Torres 
[ p7t , both to other JBD couphngs and to some particu- 
lar scalar-tensor theories with non-constant uj^cj)) chosen 
to match all current observational constraints. This al- 
lowed a study of some models which, inside the structure 
of the star, have couplings deviating greatly from the 
large value required today. The conclusion is that bo- 
son star models can exist in any scalar-tensor gravity, 
with masses which are always smaller than the general 
relativistic case (for a given central scalar field density), 
irrespective of the coupling. 

A vital point to consider is that when one finds cos- 
mological solutions in scalar-tensor theories, the gravita- 
tional coupling is normally evolving. This has important 
implications for astrophysical objects, because it means 
that the asymptotic boundary condition for the (j) field 
is in general a function of epoch. One can then ask, as 
orig inally done by Barrow in the context of black holes 
|l8| , how the structure of the astrophysical object is af- 
fected given that the asymptotic gravitational coupling 
continues to evolve after the object forms. Two possi- 
bilities exist; either the star can adjust its structure in a 
quasi-stationary manner to the asymptotic gravitational 
constant, or it might 'remember' the strength of gravity 
at the time it formed. Barrow called this latter possi- 
bility gravitational memory. In the former case, stellar 
evolution is driven entirely by gravitational effects, while 
in the latter case objects of the same mass could differ 
in other physical properties, such as their radius. Either 
possibility has fascinatingconsequences, which we have 
already explored in Ref. [19[. 

However, which of the two scenarios is correct remains 
unknown, either for black holes or boson stars. Since bo- 
son star solutions are non-singular, they appear to offer 
better prospects for determination the actual behaviour. 
Consequently, it is important at this stage to have a com- 
plete description of boson star models at different eras of 
cosmic history, which may be used later either as an ini- 
tial condition for, or to compare with the output of, a 
dynamical evolutionary code. 

Recently, two other works have been presented con- 
cerning scalar-tensor gravity effects on equilibrium boson 
stars. In the first of them. Comer and Shinkai stud- 
ied zero and higher node configurations for the Damour- 
Nordtvedt approach to scalar-tensor theories . They 
also studied stability properties of boson stars both at 
the present time and in the past. They concluded that 
no stable boson stars exist before a certain cosmic time, 
due to all possible configurations possessing a positive 
binding energy. This result appears surprising, as the 
boson stars should have no particular awareness of the 



present value of the gravitational coupling, and it would 
appear a great coincidence that the transition between 
instability and stability should occur at a recent cosmic 
epoch. In fact, their result has already been questioned 
by Whinnett in a detailed discussion of the mean- 
ing of the boson star mass in scalar-tensor theories. Our 
results also indicate that boson stars may form and be 
stable at any cosmic epoch. Finally, the dynamical for- 
mation of boson stars was analyzed in Ref. , where a 
similar behaviour to that of general relativity was found. 

In this paper, we aim to provide a comprehensive study 
of equilibrium configurations of boson stars, emphasizing 
their characteristics, such as mass and radii, at different 
moments of cosmic history. We shall also study, using 
catastrophe theory, their stability properties. As seen in 
Ref. ||l^, the features of JBD and general scalar-tensor 
boson stars do not differ much. Hence, we shall concen- 
trate only on JBD boson stars, examining the dependence 
on LU. Finally, we shall test whether the Brans-Dicke 
scalar can induce any change in the stability properties 
even for extreme values of Newton's constant. 

The organization of the rest of this work is as follows. 
In the next section we briefly introduce the formalism, 
following Ref. |17|. The following section will analyze 



some recently-proposed mass functions for JBD boson 
stars, and justify our choice for this work. We shall also 
comment on the use of catastrophe theory in the study of 
stability properties. Finally, the results of our numerical 
simulations will be given in Sec. V and our conclusions 
will be stated in Sec. VI. 



II. FORMALISM 

First we derive the equations corresponding to a gen- 
eral scalar-tensor theory. The action for these general- 
ized JBD theories is 
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Here g^i/ is the metric, R the scalar curvature, (p the 
Brans-Dicke field, and £„! the Lagrangian of the matter 
content of the system. 

We take this to be the Lagrangian density of a 
complex, massive, self-interacting scalar field ^. This 
Lagrangian reads as: 



(2) 



The U{1) symmetry leads to conservation of boson num- 
ber. Varying the action with respect to g^'^ and 4> we 
obtain the field equations: 
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where T^i, is the energy-momentum tensor for the matter 
fields and T its trace. This energy-momentum tensor is 
given by 



m^\^\^ + -m''). (5) 



Commas and semicolons are derivatives and covariant 
derivatives, respectively. The covariant derivative of this 
tensor is null. That may be proved either from the field 
equations, recalling the Bianchi identities, or by intuitive 
arguments such as the minimal coupling between the field 
(j) and the matter fields. This implies. 



mV-A|V'PV'* =0. 



(6) 



We now introduce the background metric, correspond- 
ing to a spherically-symmetric system which is the sym- 
metry we impose upon the star. Then 



-B{r)dt^ + A{r) dr^ +r^dn^ 



(7) 



We also demand a spherically-symmetric form for the 
scalar field describing the bosonic part and we adopt a 
form consistent with the static metric. 



ij^ir, t) = x(r) exp [— itui]. 



To write the equations of structure of the star, we use a 
rescaled radial coordinate, given by 

X = mr . (9) 

From now on, a prime will denote a derivative with re- 
spect to the variable x. We also define dimensionless 
quantities by 
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(10) 



— 1/2 

where mp\ = Gq is the present Planck mass. Note 
that our dimensionless variables are defined with respect 
to our observed Planck mass, regardless of whether or 
not that corresponds to the Planck mass at that time. 
Our observed gravitational coupling implies <i> = 1.^ In 
order to consider the total amount of mass of the star 
within a radius x we change the function A in the metric 
to its Schwarzschild form. 



A{x) = 1 - 



2M{x) 
X $(oo) 



(11) 



This expression defines M{x). The issue of mass defini- 
tions in JBD boson stars will be examined more deeply in 
the following section. Note that a factor appears 
in Eq. (|l^) . This is crucial to obtain the correct value of 
the mass, which from comparing this to the asymptotic 
form of the JBD-Schwarzschild solution is given by 



M,tar = M(oo) , 
m 



(12) 



for a given value of m.|^ 

With all these definitions, the equations of structure 
reduce to the following set: 
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To solve these equations numerically, we use a fourth- 
order Runge-Kutta method, for which details may be 
found in Ref. ||l^ . In general relativity, the possible equi- 
librium solutions are entirely parametrized by the central 
value of the boson field, cr(0). In JBD theory, one also 



* There is actually a post-Newtonian correction to this of 
order l/uj which we shall not concern ourselves with. 



^This corrects an error in Eq. (10) of Ref. [^. That error 
was typographical only and did not affect any computations 
in that paper. 
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needs to specify the asymptotic strength of the gravita- 
tional coupling, $00, or equivalently the value of $ at the 
centre of the star. 

The particle number, conserved due to the U{1) sym- 
metry of the ip field, is 




(17) 



where the last equality defines Noo ■ If the particles com- 
prising the star were widely separated, their mass would 
be miVstar- One can therefore define a binding energy, 
BEstar = Mstar " "i^star, and a ucccssary, though not 
sufficient, condition for the star to be stable is that the 
binding energy be negative. Considerable care is how- 
ever necessary in deciding how to define the mass which 
appears in this expression and we discuss this at 
length in the next Section. It is normally convenient to 
consider a dimensionless binding energy, defined by 



BE = Af(oo) - mN^ 



(18) 



Finally, to get a feeling for the possible rate of variation 
of we consider the solution corresponding to homoge- 

HH]: 



neous matter-dominated cosmologies, which is 
$(i) oc ^2/(4+3^) oc a^/(^+^') . 



(19) 



At the current limit w — 500, the variation in $ since 
matter-radiation equality at around Zcq = 24 000 Cloh^ is 
a couple of percent. During radiation domination $, and 
hence G, is constant. 



III. MASS DEFINITIONS 

The definition of mass in scalar-tensor theories is a 
subtle one, which has recently been examined in detail 
by Whinnett |2|]. When one leaves the security of gen- 
eral relativity, one first has to worry about which confor- 
mal frame one should work in, either the original Jordan 
frame as given in Eq. (^) , or the Einstein frame obtained 
by carrying out a conformal transformation to make the 
gravitational sector match general relativity. Addition- 
ally, while in the Einstein frame all reasonable definitions 
coincide, in the Jordan frame they do not. 

Whinnett studied three possible definitions. He found 
huge differences for w = — 1, but the three definitions 
approach each other in the large uj limit, as one expects 
since they coincide in general relativity. These are to be 
compared with the rest mass, which is just the particle 
number multiplied by the particle mass. The definitions 



where p is defined as the right-hand side of Ein- 
stein's timelike equation. This corresponds to 
the ADM mass in the Jordan frame. It is the 
commonly-used definition of mass and in the limit 
r — > cxD coincides with Mgtar defined in Eq. dl2). 



• The Keplerian mass, given by, 

mK[rj = "2" ■ 

• The Tensorial mass, given by, 

mT[r) = r" 



2<pr + r'^cj)' 



(21) 



(22) 



The last two definitions are orbital masses. A non-self- 
gravitating test particle in a circular geodesic motion in 
the geometry of Eq. (|^) moves with an angular velocity 
given by 



(23) 



as measured by an observer at infinity pq |. Then, ap- 
plying Kepler's third law, the mass of the system can be 
obtained by making, 



M(oo) = lim 



dip 
It 



(24) 



So, the Keplerian mass is Kepler's third law mass in the 
Jordan frame, whereas the Tensorial mass is Kepler's 
third law mass in the Einstein frame. In the Einstein 
frame all mass definitions coincide, so the Tensorial mass 
is also the Einstein frame ADM mass. 

These definitions differ impressively for the cu = —1 
case, and, in general, for low values of u; |2^. Then, of 
course, it becomes very important to have a correct de- 
scription of the stellar mass, because it will decide stabil- 
ity properties and binding energy behaviour. For the case 
LD = —\, the Keplerian mass would lead to positive bind- 
ing energy for all values of central density, suggesting that 
every solution is generically unstable. The Schwarzschild 
mass would instead lead to negative binding energies for 
every value of central density, suggesting that every so- 
lution is potentially stable, even for large values of cr(0).^ 
This leads one to feel that neither of these two masses 
is likely to be the correct one to use in the binding en- 
ergy calculation. Further, it is the Tensorial mass which 



The Schwarzschild mass, given by 



m(r) = 47r / pr dr , 
Jo 



(20) 



•'■In fact, for small values of central density the Schwarzschild 
mass becomes negative for low ui. This might indicate that a 
classical wormhole can form, in much in the same way as the 
solution presented in Ref. [ p7[ . 
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peaks (as a function of central density) at the same lo- 
cation as the rest mass, an important property in the 
general relativity case |^ which is crucial in permitting 
the application of catastrophe theory to analyze the sta- 
bility properties. This property presumably originates 
from the Tensorial mass being the Einstein frame ADM 
mass, though we have no mathematical proof at present. 
There is therefore a strong case towards the adop- 
tion of the Tensorial mass as the real mass of the star, 
especially for the strong field cases of low lo values. 

For the simulations we analyze in this work, we have 
computed both the Schwarzschild and the Tensorial 
mass. As expected, we find that for large w values, which 
are the ones in which we are interested, the difference is 
negligible; every graph we plot is unchanged if we replace 
the Schwarzschild mass by the Tensorial mass. Hence, for 
reasons of numerical simplicity we actually compute the 
Schwarzschild mass, as it is directly obtained from the 
set of differential field equations. 

IV. STABILITY ANALYSIS USING 
CATASTROPHE THEORY 

Catastrophe theory provides a very direct route to the 
stability properties of boson stars |29[| . The technique 
was described in a review of the stability of solitons [B0| , 
in which it was shown that the identification of conserved 
quantities of a physical system is sufficient for the deter- 
mination of stable und unstable solitons. In the case of 
boson stars, we are dealing with nontopological solitons 
which are characterized by mass and particle number, the 
only conserved quantities of this theoretical model. For 
every central value of the scalar field, there is a unique 
value for the mass and particle number. By drawing 
the conserved quantities against each other, the so-called 
bifurcation diagram is created. If cusps are present in 
this diagram, one can immediately read off the stable 
and unstable states. Starting with small central densi- 
ties where mass and particle number is also small, one 
assumes that these stars are stable (against small radial 
perturbations). If, as the central density is increased, one 
meets a cusp, the stability of the boson star changes from 
stable to unstable if the following states — the branch 
as a whole — have higher mass. This method is applied 
again at every succeeding cusp. Should it be that at some 
cusp the masses beyond the cusp are smaller, then the 
state changes from unstable to stable. The reason behind 
this method is that the cusp is a projection of a saddle 
point of a Whitney surface j2^]. The curve leading to 
the cusp consists of projections of fold points; fold points 
and cusps are the singularities of the Whitney surface, 
just the points recognizable in the bifurcation diagram. 
The fold points are the projection of maxima and min- 
ima of Whitney's surface; maxima determine unstable 
solutions while minima govern stability. 

The method of catastrophe theory has also been ap- 
plied in the context of neutron stars pi , Einstein- Yang- 
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FIG. 1. Boson star equilibrium configurations for different 
values of "I>(cxd) and self-interaction. Of the two lines for each 
case, the one with the higher peak is the particle number. 
ct(0) is in the interval (0,0.75) and there are 50 models per 
curve. 

Mills black holes , and inflationary theory . More 
recently, it has been introduced for neutron and boson 
stars in scalar-tensor theories, and in particular in JBD 
theory From these theories, one can learn how 

the properties of static solutions change if the asymptotic 
value of $ changes. In the following, we show that there 
is no stability change at all within a JBD theory as 
is changed; if a star is stable, then it is for every value 
of G. The binding energy can change its sign for stars 
which are unstable, but not for stable ones. 

V. NUMERICAL RESULTS 

First, we plot the equilibrium configuration diagrams 
for different values of the central density and asymptotic 
gravitational constant. In Fig. la (A = 0), we have 50 
models with central density in the interval (0, 0.75), with 
no self-interaction. Fig. lb shows the same, but with 
A — 100. We recognize that at fixed central density cr(0), 
the mass and particle number increase from earlier times 
($((X)) = 0.95) to later times ($(00) = 1.05). If we 
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FIG. 2. Typical curves for boson star masses as a function 
of 3>(0), with $(oo) = 0.95. Note the narrow x-axis range. 

draw the mass against the central value of the JBD field 
$(0), we find a loop, see Fig. ^. The curve starts at the 
flat spacetimc solution (<i> = constant everywhere and 
zero mass), reaches the maximum at the same value of 
central scalar field as it reached the maximum of Fig. 1 , 
cf. | pO| , makes a turn, and eventually reaches smaller 
$(0) values. Stable stars are characterized by $(0) > 
$(00), i.e. G(0) < G(oo) where G is a function of r. 
Unstable stars can have G(0) greater than or less than 
G(oo). There are two solutions for <i>(0) — $(00): first, 
the flat spacetime solution and, secondly, an unstable 
boson star. The same characteristic curve is to be found 
for different values of the asymptotic G. 

Figs. |l| and || give us, in form of (fT(0), 'i^(O)), the com- 
plete information about the initial characteristics of a 
boson star at a certain 'time', characterized by the con- 
stant <I>(oo). 

For the investigation of stability, rather than the bi- 
furcation diagram (M, TV) we use the analogous figure of 
binding energy against the particle number. Fig. |^. It 
shows us that the stars with small central densities have 
negative binding energies. In Fig. 3b (A = 100) we see 
two cusps: the first one corresponds to the maximum of 
Fig. 1 and the second to the minimum. For A = we 
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FIG. 3. The boson star binding energy as a function of 
the number of particles. The figure depict several curves for 
difi'erent <l>(oo) and A values, the coupling is a; = 400. 

did not go to high enough central densities to see the sec- 
ond cusp. The first cusp has negative binding energy and 
the other has positive binding energy. Stars with central 
densities from zero to the first cusp belong to projections 
of minima within a Whitney surface, i.e. they are stable. 
Beyond the cusp, one radial perturbation mode is be- 
coming unstable, and at the second cusp a second mode 
becomes unstable. 

To study the influence of the changing gravitational 
coupling on the exact position of the cusp we made a 
high resolution study around the position of the cusp for 
the <I>((X)) = 1 model. As can be seen from from Fig. 1, it 
is at (t(0) ~ 0.27. We then did simulations in the interval 
cr(0) e [0.265,0.275], with eleven models in that range, 
shown in Fig. 4. For $(00) = 1.05 the mass and the 
number of particles as a function of cr(0) are increasing 
functions. That means all the models are in the stable 
branch. For $(00) = 0.95 instead, mass and number of 
particles are decreasing functions, locating all models in 
the unstable branch. In the case of our present gravi- 
tational coupling, the cusp appears within the interval. 
Thus, going from the future to the past {${00) — 1.05 to 
$(00) — 0.95) models with a given central density move 



6 




o 0.266 0.268 0.27 0.272 0.274 



FIG. 4. Boson star masses for a reduced interval of values 
of (j(0) around the cusp at three different 'times' <l>(oo) — 1.05 
(top), $(oo) = 1.0 (middle), and ${oo) = 0.95 (bottom) 
[uj = 400]. 

from the stable to the unstable branch. This agrees with 
the analytic prediction that in the general relativity limit 
one should find cr^ax ^{oo). The movement of the cusp 
is much the same as Comer and Shinkai reported in [ pO[ , 
except for one important point. They found no cusp at 
all for times well before the present, meaning that they 
did not find any stable star in the past. On the contrary, 
we have found that the cusp moves backwards in cr(0), 
but it is still there, see Fig. 3. We believe that their con- 
clusion derives from the use of a wrong mass definition. 

In addition, we have calculated solutions with constant 
central scalar field values at different 'times' <&(oo), also 
taking into account very small values of <I>(oo) which are 
unphysical, see Fig. |^. This figure represents a bifurca- 
tion diagram with respect to $(cxd). It is evident that no 
cusp is present, so no stability change occurs. 

Because a boson star has no clearly-defined surface, 
but rather an infinite exponentially-decreasing atmo- 
sphere, several radius definitions are in use. We apply 
here the common one, the radius which encloses 95% of 
the total mass. Fig. ^represents the mass against the ra- 
dius. The diagram shows that solutions with small cen- 
tral densities have large radii. Then, with growing central 
densities, the mass increases while the radius decreases. 
The maximum in this diagram is the most centrally-dense 
stable star solution. The radius of the maximal mass bo- 
son star remains roughly the same, but the mass corre- 
sponding to a given central density grows with time, pro- 
ducing a denser star. The increase of the self-interaction 
constant A gives larger radii as one expects from a rcpul- 




0.0 I ' ' ' ' 1 

0.5 1.0 1.5 2.0 2.5 

Phi(infty) 



1.5 - 




0.0 I 1 

0.5 1.0 1.5 

N 

FIG. 5. The behaviour of masses and particle numbers for 
extreme values of the gravitational asymptotic constant. The 
model taken has A = 0, cr(0) =0.1 and ui = 400. 

sive force. Compare this with similar results for neutron 
stars (Figure 7 in Ref. |^^) and for general relativistic 
boson stars (Figure 3 in Ref. |3l|). 

So far, we have recognized that the boson stars are 
denser the larger $(oo) is. The reason can be understood 
as a deeper gravitational potential, expressed by an in- 
crease in the difference of $(0) and <&(oo), see Fig. |^. 

For a fixed value of cr(0), the behaviour of the bind- 
ing energy, the radius, and 6^ (the difference between 
the central and the asymptotic value of the Brans-Dicke 
scalar) are all plotted in Fig. ^. 

In Fig. ^ we show the dependence of equilibrium con- 
figurations on uj. To do this we plot the binding energy 
behaviour in the interval cr(0) G (0,0.3) for two values 
in the asymptotic effective gravitational constant. The 
value of oj is in the range (50, 50000). The upper curves 
in both diagrams correspond to ui — 10000 and 50000 
and match each other exactly. The upper panel shows 
models with <I>(c») = 0.98, while the other has our ob- 
served gravitational strength. This shows that, when u; 
tends to infinity, a general relativity like solution — with 
a different value for Newton's constant — is obtained. 
Recall that even with the strong limit on to valid today, 
we could have an evolving uj{(f>) which is much smaller in 
the past, and so small values of the coupling parameter 
may also be meaningful. 
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FIG. 6. The radius of equilibrium boson star configura- 
tions for different values of the self-interaction and central 
density in the range (0,0.75). 

VI. CONCLUDING REMARKS 




ct(O) 

FIG. 7. The behaviour of the difference between the cen- 
tral and the asymptotic value of the Brans-Dicke scalar as a 
function of ct(0) for different values of the effective gravita- 
tional constant. Note the highly expanded j/-axis. 



obtained in this paper can be regarded as of a general 
nature. Astrophysics should be unambiguously sensitive 
to the underlying theory of gravity, especially on cosmo- 
logical times scales. It is in this framework, perhaps, 
where a crucial test of gravity could arise. 
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In this paper, we have thoroughly analyzed static bo- 
son star configurations in the framework of the Jordan- 
Brans Dickc theory of gravitation. Wc studied their 
equilibrium and stability properties in the present as well 
as for other cosmic times, in the past or in the future. 
Stable boson stars may exist at any epoch, with stabil- 
ity depending on the value of central density. Together 
with this, a number of new physical features have been 
displayed concerning the radius-mass relation, the be- 
haviour of the difference between the central and asymp- 
totic value of $, the dependence on the structure upon 
the coupling parameter and other properties. This con- 
figurations can be used either to compare with the output 
of a numerical evolution code, or as the input into one. 
We expect that such a study will shed light on which 
scenario of the gravitational memory phenomenon might 
occur in practice. Whichever it might be, it is very likely 
that the same phenomena could also occur for fermionic 
stars, such as white dwarfs. In this sense, the results 
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